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ABSTRACT 

We study the halo mass accretion history (MAH) and its correlation with the internal 
structural properties in coupled dark energy cosmologies (cDE) . To accurately predict 
all the non-linear effects caused by dark interactions, we use the COupled Dark Energy 
Cosmological Simulations (CoDECS). We measure the halo concentration at z = and 
the number of substructures above a mass resolution threshold for each halo. Tracing 
the halo merging history trees back in time, following the mass of the main halo, we 
develope a MAH model that accurately reproduces the halo growth in term of M 2 oo in 
the ACDM Universe; we then compare the MAH in different cosmological scenarios. 
For cDE models with a weak constant coupling, our MAH model can reproduce the 
simulation results, within 10% of accuracy, by suitably rescaling the normalization of 
the linear matter power spectrum at z = 0, eg- However, this is not the case for more 
complex scenarios, like the "bouncing" cDE model, for which the numerical analysis 
shows a rapid growth of haloes at high redshifts, that cannot be reproduced by simply 
rescaling the value of <r$. Moreover, at fixed value of as, cold dark matter (CDM) 
haloes in these cDE scenarios tend to be more concentrated and have a larger amount 
of substructures with respect to ACDM predictions. Finally, we present an accurate 
model that relates the halo concentration to the time at which it assembles half or 4% 
of its mass. Combining this with our MAH model, we show how halo concentrations 
change while varying only er s in a ACDM Universe, at fixed halo mass. 
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1 INTRODUCTION 

Understanding the formation and evolution of structures in 
the Universe is one of the main goals of present cosmologi- 
cal studies. Following the standard scenario, the formation 
of cosmic structures up to protogalactic scale is due to the 



gravitational instability of the dark matter (DM) (Frenk 



et al.|1983||Davis et al.|1985l|White|1988||Frenk et al.|1990 



Springel et al.||2005| 2008 1 . When a density fluctuation ex- 
ceeds a certain value, it collapses forming a so-called DM 
halo. The small systems collapse first in a denser Universe 
and then merge together forming the larger ones, thereby 
giving rise to a hierarchical process of structure formation. 
In this scenario, galaxy clusters are located at the top of 
the merger history pyramid and today represent the largest 
virialized objects in the Universe. The formation of luminous 
objects happens when baryons, feeling the potential wells of 
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the DM haloes, fall inside them shocking, cooling and even- 
tually forming stars (White & Rccs 19781. New supplies of 
gas and galaxy mergers tend to modify the dynamical and 
morphological properties of the forming systems, and are 
closely linked to the mass accretion histories of the haloes 
they inhabit. A detailed understanding of how this mass ac- 
cretion occurs and how individual halo properties depend 
on their merger histories is of fundamental importance for 
predicting galaxy properties within the CDM theory and, 
similarly, for using observed galaxy properties (as e.g. rota- 
tion curves) to test the paradigm. 

Different definitions have been adopted in the litera- 
ture to study the halo growth along the cosmic time and 



its correlation with the halo clustering on large scales ( Gao 
|et al.||2004| |Gao fc Whlte||2007| ). From different analyses it 
has emerged that while the redshit, Zf, at which the main 
halo progenitor assembles a fraction, /, of its mass correlates 
with its global structural properties (as e.g. its concentra- 
tion, spin parameter, subhalo population, etc.), the redshift, 
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z c , at which the main halo progenitor assembles a constant 
central mass, M c , mainly correlates with the typical forma- 
tion time of stars in a halo. Considering the haloes at z = 
in the Millennium Simulation, Li et al. ( 2008[ ) have shown 
that, while zj decreases with the halo mass, z c grows with 
it, in agreement with the fact that older stellar populations 
tend to reside in more massive systems. 

Many studies conducted on standard ACDM simula- 
tions have also underlined how structural halo properties, 
like concentration and subhalo population, are related to 
the MAH ( |van den Bosch|2002| [Gao et al.|2004| |2008[ >. Less 
massive haloes, typically assembling a given fraction of their 
mass earlier, tend to host less substructures than the more 



massive ones. In particular, Giocoli et al. ( 2008 ) found that 



at fixed halo mass more concentrated haloes tend to possess 
less substructures because forming at higher redshift than 
the less concentrated ones. Extending these results to the 
framework of the assembly bias we would expect - at fixed 
halo mass - more concentrated, more relaxed, and less sub- 
structured haloes to be on average more biased with respect 
to the DM. 

Analytical models of the DM density distribution, based 
on the halo mod el ( |Seljak|2000| |Cooray fc Sheth|2002| |Gio~ 



coli et al. 2010), require knowledge of the halo mass func- 



tion, density profile, concentration and subhalo population; 
as does the halo occupation distribution (HOD) approach 
to describe the galaxy and the quasar luminosity function 
and their bias (|Moster et al.||20lo| |Shen et alT||2010| |Cac-| 
ciato et al.||2012 i. At the same time, these results are use- 
ful to model the weak and strong lensing signals by large 
scale structures (jBartelmann fc Schneider|2001[ ) and clusters 



I Giocoli et al.|2012a|) in standard and non-standard cosmolo- 



gies. Since many cosmological models have been proposed as 
a possible alternative to the standard concordance ACDM 
scenario, it is natural to investigate whether these results 
can also be extended to non-standard cosmologies. In this 
paper, we will focus on a range of DE scenarios character- 
ized by a non-vanishing coupling between the DE field and 



CDM particles, which go under the name of cDE ( Wetterich 
1995{ |Amendola|2000| [2004 1 



This paper is divided in two parts. In the first part, 
we investigate the impact of a number of cDE models on 
the formation and accretion histories of CDM haloes (S|2] - 
Sj5|. To tackle this point, we make use of two theoretical 
techniques, one analytic and one numeric. The former one 
allows us to predict the MAH through a generalized version 



of the Press & Schechter ( 1974 ) formalism. The numerical 



approach aimed at properly describing all the non-linear ef- 
fects at work. Here, we consider the mock halo catalogues 
extracted from the CoDECS simulations and compare them 
to our analytic MAH predictions. In the second part of the 
paper (|j6j, we exploit both the above theoretical methods 
to predict the internal structural properties of CDM haloes 
in these cosmological scenarios. In particular, we derive the 
halo concentration-mass relation, that is shown to provide 
a key observable to discriminate between cDE cosmologies 
and, in some cases, to remove the degeneracy with the nor- 
maliazion of the power spectrum. 

After a general introduction to the cDE cosmologies 
analyzed in this work in ^2j we describe our numerical N- 
body simulations in S|3] Details on the method we use to 
climb the halo merging history trees are given in Q while 



our theoretical predictions, both analytical and numerical, 
on the halo MAH are provided in [J5] In !|6] we investigate 
the halo internal structural properties. Finally, in SjTJwe draw 
our conclusions. 



2 THE CDE MODELS 



cDE models have been proposed as a possible alternative 
to the cosmological constant A and to standard uncou- 
pled Quintessence models ( Wetterich|1988 Ratra & Peebles 
1988| ) to explain the observed accelerated expansion of the 
Universe (|Riess et al.|1998||Perlmutter et al.|1999||Schmidt 



et al. |1998[ ), as they alleviate some of the fine-tuning prob- 
lems that characterize the latter scenarios. In order to be 
viable, cDE models need to assume a negligible interaction 
of the DE field to baryonic particles ( jDamour et al.|[l990 l, 
since the long-range fifth-force that arises between coupled 
particle pairs as a consequence of the interaction with the 
DE field would otherwise violate solar-system constraints 
on scalar-tensor theories (see e.g. Bertotti et al.||2~003 Will 
|2005[ > . Consequently, a large number of cDE models charac- 
terized by a direct interaction between DE and CDM (see 



e.g.|Wetterich|1995||Amendola|2000l[200^[FarFar fc Peebles 



2004[|Caldera-Cabral et al.|2009||Koyama et al.|2009||Baldi 



2011b I or massive neutrinos (see e.g. Wetterich|2007 Amen- 



dola et al. 2008) have been proposed in recent years. The 



basic properties of cDE models and their impact on struc- 
ture formation have been thoroughly discussed in the lit- 
erature. For a self-consistent introduction on cDE scenarios 



we suggest for example the recent reviews Tsujikawa (|2010 1 



|De Felice fc Tsujikawa| ( |2010[ ) ; |Amendola et al.| ( |2012| - Sec- 
tion 1.4.3). For the aims of the present work, it is sufficient to 
summarize the main features that characterize cDE models 
in general, and the specific realizations of the cDE scenario 
that will be investigated here. 

In general, cDE models are characterized by two free 
functions that fully determine the background evolution of 
the Universe and the linear and non-linear growth of density 
perturbations. These are the self-interaction potential V(4>) 
and the coupling function P{tf>), where is a classical scalar 
field playing the role of the cosmic DE. The coupling func- 
tion /3((f>) determines the strength of the interaction between 
DE and CDM (in the case investigated here) and directly 
affects the evolution of density perturbations through: i) a 
long-range attractive fifth-force of order /3 2 times the gravi- 
tational strength, and ii) a velocity-dependent accele ration 



on coupled CDM particles proportional to <f>j3 (see e.g. Baldi 
|2011a| ). Differently from the fifth- force term, which is always 
attractive regardless of the sign of the coupling function and 
of the dynamical evolution of the DE scalar field, the friction 
term can take both positive and negative signs depending on 
the relative signs of the scalar field velocity <f> and the cou- 
pling function /3((/>). Such feature can have very significant 
consequences on the evolution of structure formation in case 
the friction term changes sign during the cosmic evolution 



I Baldi 2012a Tarrant et al. 2012 1 as for the case of the 
"bouncing" cDE model (Baldi 20 12a| l investigated in this 
work. 
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Table 1. The cDE models of the CoDECS suite considered in the present work. All models have the same normalization of scalar 
perturbations at .zcmb ~ 1100 leading to a different value of erg. Besides trg and the value of the DE equation of state all the models 
share the same WMAP7 | |Komatsu et al.|2011| cosmological parameters at the present time. 

Model Potential a /3 /3i Scalar field Potential 10^(2 = 0) A s {zcmb) ct s (z = 0) 

normalization normalization 



ACDM 


v(fl 


= A 












A = 


0.0219 


-1.0 


2.42 X 10~ 9 


0.809 


EXP003 


v(4>) 


= Ae-™^ 


0.08 


0.15 





4>(z = 0) = 





A = 


0.0218 


-0.992 


2.42 X 10" 9 


0.967 


EXP008e3 


V(4>) 


= Ae- a ^ 


0.08 


0.4 


3 


4>{z = 0) = 





A = 


0.0217 


-0.982 


2.42 X 10" 9 


0.895 


SUGRA003 


V(4>) 


= A0- Q e<* 2 / 2 


2.15 


-0.15 





(j>(z — > 00) = 




A = 


0.0202 


-0.901 


2.42 X 10~ 9 


0.806 



3 NUMERICAL SIMULATIONS 

We make use of the publi c halo and su bhalo catalogues of 
the CoDEC^j] simulations ( |Baldi||2012b| - the largest suite 
of cosmological N-body simulations of cDE models to date 
- to follow the accretion histories of CDM haloes by build- 
ing the full merger trees of all the structures identified at 
z — up to z — 60. In our analysis, we will consider the dif- 
ferent combinations of the potential and coupling functions 
that are available within the CoDECS suite of cDE scenarios, 
defined by the following general expressions: 



• Exponential potential ( Lucchin fc Matarrese|1985 Wet- 
|terich|1988| : 

V{cj>) = Ae- a +; (1) 

• SUGRA potential ( |Brax fc Martin|1999[ ): 

V(4>) = A^ a e^ /2 ; (2) 

• Constant coupling (Wcttcrich 1995; Amcndola 2000): 

P{4>) =p = const. ; (3) 
Time-dependent coupling (Amcndolaj |2004| |Baldi| 



2011b) 



0(4>) = Poe 01 < 



(4) 



In particular, while standard cDE models are characterized 
by "runaway" potentials (like e.g. the exponential poten- 
tial of Eq. [T]) and a constant coupling, time-dependent cDE 
models feature the same type of potentials with a coupling 
function that evolves with the scalar field itself. Finally, the 
combination of a "confining" potential like the SUGRA po- 
tential of Eq. [2] with a negative constant coupling character- 
izes the so-called "bouncing" cDE scenario (Baldi 2012a). 



A summary of the cDE models that are investigated in the 
present work, with the corresponding parameters, is given in 
Table [l] For such scenarios, we will make use of the publicly- 
available data of the L-CoDECS series to derive halo accre- 
tion histories in the different cosmologies. 

The L-CoDECS runs are collisionless N-body simula- 
tions carried out on a periodic cosmological box of 1 Gpc//i 
aside, filled with 2 x 1024 3 equally sampling the coupled 
CDM and the uncoupled baryon fluids. The baryons are 
treated as a separate family of collisionless particles, and 
no hydrodynamic treatment is included in the simulations. 
The mass resolution at z = is m c = 5.84 x 10 10 M e /h 
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and m b = 1.17 x 10 10 M Q /h for CDM and baryon parti- 
cles, respectively, while the gravitational softening is e g — 20 
kpc//i. All the models virtually share the same initial con- 
ditions at the redshift of last scattering z\ a ~ 1100 and are 
consequently characterized by a different value of erg, due to 
their different growth histories (as summarized in Table [T]). 
Except for the different values of erg and of the DE equation 
of state parameter w<j,, all the models share the same cos- 
mological parameters at z = 0, consistent with the WMAP7 
results (Komatsu et al. (2011 1, see Table|2|. This provides a 
self-consistent set of cosmologies that can be directly com- 
pared with each other and tested with present and future 
observations. The viability of these models in terms of CMB 
observables has yet to be properly investigated, in particu- 
lar for what concerns the impact of variable-coupling and 
"bouncing" cDE models on the large-scale power of CMB 
anisotropics . Although such analysis might possibly lead to 
tighter bounds on the coupling and on the potential func- 
tions than the ones allowed in the present work, here we are 
mainly interested in qualitatively understanding the impact 
of cDE scenarios on the formation history of CDM haloes at 
late times, and we deliberately choose quite extreme values 
of the models parameters in order to maximize the effects 
under investigation. 

The CoDECS simulations have already been used for 



several investigations. In particular, Lee & Baldi (20111 ex- 



ploited these data to study the pairwise infall velocity of 
colliding galaxy clusters in cDE models, demonstrating that 
DE interactions can significantly enhance the probability 
of high- velocity collisions. |Beynon et al. | ( |2012[ ) provided 
forecasts for the weak lensing constraining power of future 



galaxy surveys, while Cui et al. (20121 exploited the same 
data to test the universality of the halo mass function. Fi- 
nally, the CoDECS data have been used to study how DE 
interactions modify the halo clustering, bias and redshift- 
space distortions ( jMarulli et al.||2012|), and how they can 
shift the baryonic acoustic oscillations JVera Cervantes et al. 



2012) 



4 THE HALO MERGER HISTORY TREE 

For each simulation snapshot, we identify haloes on the fly 
by means of a Fried-of-Friend (FoF) algorithm with link- 
ing parameter b = 0.2 times the mean interparticle separa- 
tion. Within each FoF group we also identify gravitationally 



bound substructures using the SUBFIND algorithm ( Springel 
|et al.|2001b"t . SUBFIND searches for overdense regions within 
a FoF group using a local SPH density estimate, identifying 
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Table 2. The set of cosmological parameters assumed for all the 
models included in the CoDECS Project, consistent with the 7 
year results of the WMAP collaboration for CMB data alone ( |Ko-| 
|matsu et al.|2011) . 



Parameter 



Value 



H 

fiCDM 



70.3 km s _1 Mpc" 
0.226 
0.729 
0.809 
0.0451 
0.966 



substructure candidates as regions bounded by an isoden- 
sity surface that crosses a saddle point of the density field, 
and testing that these possible substructures are physically 
bound with an iterative unbinding procedure. For both FoF 
and SUBFIND catalogues we select and store systems with 
more than 20 particles, and define their centres as the posi- 
tion of the particle with the minimum gravitational poten- 
tial. It is worth to notice that while the subhaloes have a 
well defined mass that is the sum of the mass of all particles 
belonging to them, different mass definitions are associated 



to the FoF groups (Lukic et al. 20091. We define as Mf f 



the sum of the masses of all particles belonging to the FoF 
group, M200 the mass around the FoF center enclosing a 
density that is 200 times the critical one, and M„; r as the 
one enclosing the virial overdensity A„j r , as predicted by the 
spherical collapse model. We notice that in the ACDM Uni- 
verse at z = generally it is true that M200 < M V i r < Mf f 



( Eke et al.ft996| |Bryan fc Norman[ 1998 ) . 

It is interesting to notice that the definition of M„; r is 
cosmology dependent, since it is related to the A v i r defi- 
nition derived from the spherical collapse model. Generally 
this quantity is not directly obtainable for any arbitrary DM 



and dark energy model. For example, Pace et al. (2010) have 



shown how to compute A„i r from the non-linear differen- 
tial equation for the evolution of the density contrast for 
dynamical and early DE cosmologies. For ACDM models 
many fitting functions are available that depend on the red- 



shift evolution of the DM content in the Universe (Navarro 
et al. 1997] Bryan & Norman 19981, while no simple formu- 



lae exist A„j r within the cDE models investigated in this 
work. Since M200 is cosmology independent, we will adopt 
this mass definition in what follows. 

For each subhalo, starting from redshift 2 = 0, we fol- 
low back its merger history tree by requiring to have at least 



one descended subhalo at the previous snapshot (Boylan 



Kolchin et al. 20091. In order to trace back in time the 



growth in terms of M200, we consider only the dominant 
subhalo of each FoF group to which is associated the value 
of A/200 of the group. In this process we link together among 
the different simulation snapshots the main progenitor, de- 
fined as the dominant subhalo that donate the largest num- 
ber of particles between two consecutive snapshots and the 
satellites that represent the progenitor haloes accreted onto 
the main progenitor branch of the tree ( Tormen et al.|2004 |. 



5 THE MASS ACCRETION HISTORY 

Before introduce the results of the other cDE models, in this 
section, we review the model of the mass accretion history 



proposed by Giocoli et al. (2012b I for a ACDM model and 



modify it to be consistent with the M200 halo mass defini- 
tion. This simple model allows us to derive the generalized 
formation redshift distribution, defined as the redshift at 
which the main halo progenitor assembles a fraction / of 
its M200 mass at z = 0, or more in general at any consid- 
ered redshift zq. Then, from such distribution we can com- 
pute the halo mass accretion history. Since haloes form as a 
consequence of gravitational instability processes that start 
within the initial DM density field and then grow together 
via merging events, a few important consequences are: (i) 
the MAH depends on the initial density fluctuation field - 
and so on the initial matter power spectrum; (ii) it also de- 
pends on cosmology through the background expansion of 
the Universe, (iii) at a given redshift more massive haloes 
grow faster - since they sit on the highest peaks of the den- 
sity fluctuation field; (iv) at a given mass, the higher the 
redshift is the higher the halo mass growth rate is. These 
points will become more clear along the discussion of the 
results presented in the following sections. 



5.1 Mass Accretion History in a ACDM Universe: 



revisiting the Giocoli et al. 

M200 



(2012b) model for 



Giocoli et al. (2012b) developed a simple and accurate model 



to describe the growth of CDM haloes along cosmic time. 
Their study was conducted analyzing the merger trees ex- 
tracted from a cosmological N-body simulation of a ACDM 



Universe, the GIF2 simulation (Gao et al. 2004 Giocoli et al 



2008 2010a) 



Before describing the model, let us first introduce 
some universal quantities that will be used throughout this 
work. We define S(M) as the variance of linear fluctua- 
tion field when smoothed with a top-hat filter on a scale 
R = (3M/4ttp) 1/3 : 



S(M) 



1 

2^2 



W 2 (kR)P lin (k)k 2 dk, 



(•5) 



where p is the comoving density of the background, P{k)u„ 
the linear matter power spectrum and W(kR) the top-hat 
window function. We also define S c (z) as the initial overden- 
sity required for spherical collapse at redshift z: 



5 c {z) 



Sc.linjz) 



(6) 



where 5 ,Un(z) is the linear overdensity at redshift z and 
D + (z) the growth factor normalized to unity at the present 
time. For a CDM Universe, [Nakamura & Suto ( 1997 1 have 
presented a useful fitting function for the linear overdensity 
as a function of the matter density parameter, Q m , that can 
be read as: 

S c ,iin{z) = 1.686[1 + 0.123 log(n m (*)] (7) 

and that will be used throughout this paper. In Figure [I] 
we show the mass variance for the four cosmological models 
considered in this work. We notice that all of them have the 
same shape but a different normalization. This is due to the 
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different values of as defining the variance of linear density 
fluctuations on a scale of 8 Mpc/h: 



\{M) = 



2tt 2 



W 2 (8k)P Hn (k)k 2 dk- 



(8) 



the higher is as the higher the normalization of the mass 
variance and viceversa. 



Following the formalism by Lacey & Cole (19931, we 
can define the formation redshift Zf of a halo of mass M 
(at redshift zo) as the redshift at which the main halo pro- 
genitor assembles for the first time half of its mass. Their 
proposed redshift distribution at half-mass, in terms of uni- 
versal variables, can be read as: 



p(u>f) = 2wf erfc 



Wf 

V2 



(9) 



where Wf = (8 C (z ) - 5 c (z f ))/y/S(M/2) - S(M). It is in- 
teresting to notice that written in this way the formation 
redshift distribution is both independent of the halo mass 
M and of th e redshift zo. In a more general way |Nusser fc| 



Sheth ( 1999 1 have developed a formalism to describe the 
redshift at which the main halo progenitor assembles a frac- 
tion 1/2 / < 1 of its mass, that is: 



p(u>f) = 2wf 



cxp 



(10) 

which for / = 1/2 gives back the above expression Q. 
However, comparing these distributions with measurements 



performed on numerical simulations, Giocoli et al. (20071 
have shown that they tend to underestimate the formation 
redshift distribution. To reconcile theory and simulations 
they proposed to modify Wf — > \/0.707u>/ in the context of 
the ellipsoidal colla p se model pheth fc Tormen|1999| |2002 



|Despali et al.||2012[ ). |Giocoli et al.| ( |2012b| ) have shown that 
in numerical simulations the modified equation ( 10 1 works 
quite well for smaller values of /. They also provided a new 
function to describe the formation redshift distribution for 
any fraction < / < 1: 



p(w f ) 



afwe 



u2 /2 + a/ _ 1] 



(11) 



where a / represents a free parameter. By integrating Eq. [TT] 
one gets the following cumulative generalized formation red- 
shift distribution: 



wi/2 , 

e f + ctf 



(12) 



Since the previous equation can be inverted, it is possible 
to analytically compute the median value Wf = Wf given 
by definition when P(> Wf) = 1/2; this gives the median 
redshift zj = 8~ 1 (zf) at which the main halo progenitor 
assembles a fraction / of its mass M : 



where 



S c (zf) = 5 c {z ) + w f VS(fM) - S(M) , (13) 



Wf 



y/2 In (a s + 1) 



(14) 



In Figure [2] we show the cumulative generalized for- 
mation redshift distribution in terms of the rescaled vari- 
able Wf when the main halo progenitors assemble 90%, 
50%, 10% and 1% of their M 200 mass at z = 0. The data 
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Figure 1. The mass variance 8 (see equation ^ as a function 
of the halo mass. The different curves show the results for the 
different linear power spectra of the models considered in this 
work. Their shapes are independent of the cosmological model 
considered while the normalizations change depending on the lin- 
ear power spectrum normalization and so on the corresponding 
as. 



points represent the measurements in the ACDM simula- 
tion considering all haloes at zo = with mass larger than 
M200 > 5 x 10 12 M Q //i (which means that systems are re- 
solved with at least 40 particles) and that never exceed along 
their growth histories more than 10% of their mass at the 
present time. This ensures to consider only haloes that grow 
mainly hierarchically and do not fragment as a consequence 
of violent merging events. With respect to |Giocoli et al.| 
| |2012b| >, in this case haloes are followed back in time con- 
sidering M200 as mass definition instead of M„; r The solid 



curves show the best-fit model of Eq. ( 12 1, while the dashed 
curves represent equation ( | 10[ > modified as suggested by |Gio-| 
coli et al.l (120071). 



From the figure we notice that - since at z = the value 
of M200 is smaller than that of M V i r - following the merger 
tree back in time in terms of M200 results in a value of Zf 
typically higher than that obtained by following back the 



tree in terms of M v i r (see Giocoli et al. 2012b or the Ap- 



pendix [X] where we show the same cumulative distribution 
of formation redshift following the haloes in term of their 
virial mass for the ACDM run). 

In Figure [3] we show the best fit a/ as a function of 
/. The open circles show the best fit to the ACDM simula- 
tion measurement for different values of the assembled mass 
fraction, while the solid line refers to the best fit to the data 
given by: 



1.365 _ 2f 3 
a f = Joss e 



(15) 



Finally, the shaded regions represent the la and 2a contours 
of A X 2 , where X 2 («/) = EJ« - P (> ™i,f,<*f?l 



6 Giocoli et al. 2012 




Figure 2. Cumulative generalized formation distribution, as a function of wj, when the main halo progenitor assembles 90% (top-left 
panel ), 50% (top-right panel), 10% (bottom- left panel) and 1% (bottom-right panel) of its M200 mass at redshift zero. The data points 
show the measurements for the ACDM simulation considering all haloes with at least 40 particles. The dashed curve on the top panels 
refers to the model by [Nusscr & Shcth] | |1999| l - valid only for 1/2 < / ^ 1, while the solid curve represents eq. |l2l l with the best fit 
parameter aj. We recall that we consider only haloes that along their growth never exceed more than 10% their mass at the present 
time. 



5.2 Generalized Formation Redshift and Mass 
Accretion History 

The hierarchical model predicts that small haloes tend to 
form at higher redshift and then merge together forming the 
larger ones. The halo collapse happens when a density fluc- 
tuation exceeds the critical value predicted by the spherical 
collapse model. Typically, within a ACDM Universe fluctua- 
tions in a density field with a larger amplitude collapse ear- 
lier. This can be rephrased considering two identical CDM 
initial density fields but with different normalization param- 
eters erg: haloes in the Universe with higher erg will collapse 
earlier than those in the Universe with lower as- This also 
results in the fact of same mass haloes being more concen- 



trated in a Universe with a larger value of erg ( Maccid et al 



2008 1 . In what follows we will try to understand if this holds 
also for different cosmological scenarios, i.e. if the mass ac- 
cretion history measured in numerical simulations of cDE 



can be reproduced by the MAH model built for ACDM only 
by suitably changing the linear power spectrum normaliza- 
tion erg only. 

In the top panels of Figure [4] we show the median for- 
mation redshift as a function of the halo mass, considering 
/ = 0.9, 0.5, 0.1, 0.01 as assembled mass fractions. The open 
circles show the median of the measurements in the ACDM 
simulation and the shaded region encloses the first and the 
third quartiles of the distribution at fixed halo mass. On av- 
erage a present-day cluster size halo with mass of 1O 14 M0 jh 
assembles 90% of its mass at z = 0.2, 50% at 0.7 and 1% at 
approximately z = 4. The data points in the bottom pan- 
els represent the differences of the median measured in the 
three cDE models with respect to the ACDM one. In each 
panel, the three curves show the model of the Zf — M200 re- 
lation rescaled with respect to the ACDM model computed 
using the formalism described in the previous section. 
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0.4 0.6 
fraction 

Figure 3. Dependence of ctf, the free parameter in equation | |12[ ), 
on the mass fraction / = MaooCzf )/-M20o(0) assumed at the for- 
mation. The shaded regions represent the lcr and 2a contours. 
The blue curve represent the best fit relation. 



It is important to underline that in order to build the 
MAH model in the three cDE cosmologies we need to know 
the initial density fluctuation field to compute S(M), and 
the redshift evolution and the growth of the perturbations 
given by S c (z) and A+(z). For the definition of S(M), we use 
for the three cDE models their corresponding linear power 
spectra, that are obtained from the ACDM one by renor- 
malizing the latter adopting the different values of as char- 
acterizing each cDE model (see Table [TJ. Since the purpose 
of this work is to understand if the MAH of the cDE sim- 
ulations can be obtained from the ACDM by rescaling as 
only, we adopt for the definitions of both S c (z) and D+ the 
ones obtained for the latter model. From Fig. [4] we notice 
that only for / = 0.5 the three cDE results are quite in 
agreement with the MAH model where the power spectrum 
normalization changes. Among them, since the SUGRA003 
model has the same ag of ACDM, we would expect to find a 
formation-redshift relation quite similar to the latter. How- 
ever, this is clearly not the case due to the markedly different 
growth of SUGRA003 as compared to ACDM. For / = 0.9 
and / = 0.5 haloes in the SUGRA003 model form typically 
at the same redshift as in the ACDM run, while for / = 0.1 
and 0.01 the difference appears of the order of 10%. We un- 
derline also that for the two models EXP003 and EXP008e3 
the halo formation redshifts are higher for any / than those 
measured in ACDM for the same mass haloes at z = 0. 

To better understand how on average haloes assemble 
as a function of redshift, we show in Figure [5] the median 
mass accretion history for three different mass bins at z — 
by defining $ as the logarithm of the final assembled mass 
fraction: 



*(z) = log 



M 200 (z) 



M 200 (0) 



In the top panels we show the measurements for three dif- 
ferent final halo mass bins at the present time in the ACDM 
simulation. The shaded region encloses the first and the 
third quartiles of the distribution at fixed halo mass frac- 
tion. The solid curve represents the model built using equa- 



tion ( 13 \ for different values of the assembled halo mass 



fraction < / < 1. 

The model predicts for the ACDM run a MAH that is 
in agreement within a few percent with the numerical simu- 
lation measurements both down to very small values of the 
assembled mass fraction and up to high redshifts. In the bot- 
tom panels we show the percentage difference of the MAH 
measured in EXP003, EXP008e3 and SUGRA003 with re- 
spect to the one in ACDM. In the figure, the different data 
points are the same as in Figure [4] The three solid curves 
represent the difference of MAH models changing the lin- 
ear power spectrum definition with respect to the ACDM 
model. As expected, since the SUGRA003 model has al- 
most the same linear power spectrum normalization of the 
ACDM, the MAH model predicts the same growth of the 
latter, in contrast with what is measured in the numeri- 
cal simulation. For EXP003 the halo growth is quite well 
captured changing the linear power spectrum normalization 
in our model - and is therefore fully degenerate with ag. 
However, for EXP008e3 the model fails for large haloes and 
high redshift. This can be better seen in Figure [6] where we 
show the residuals of the MAH measured in the three cDE 
simulations with respect to the MAH model obtained from 



equation ( 13 \ changing the power spectrum normalization 



accordingly. The vertical dashed line in each panel shows 
the redshift at which the mass accretion history model falls 
below 5 x 10 12 Af s //i, mass above which haloes are resolved 
in the simulations with at least 40 particles. 



6 THE HALO STRUCTURAL PROPERTIES 
6.1 Halo concentrations 

One of the most important result obtained through N-body 
simulations of structure formation is that the CDM density 
distribution in collapsed haloes tends to follow a universal 
profile. Both for small and large mass haloes the density 
profile is well described by the Navarro-Frenk- White (NFW, 



Navarro et al. 19971 relation that reads as: 



p(x) 



Pa 



x(l + x) 2 



(17) 



where x = r/r s , r s = -R200/C200 is the scale radius where the 
logarithmic slope of the density profile approaches to —2, p s 
is the density enclosed within the scale radius, and C200 is 
the halo concentration parameter. Denoting with M200 the 
mass enclosed within -R200 we can write: 



M200 
47rr§ 



ln(l + C200) 



C200 
1 + C200 



(18) 



= log(/) • 



(16) 



To characterize the halo concentration in the numerical 
simulations we use the same approach adopted by |Springel| 
|et al.| ( |2008[ ) and |Cui et al.| ( |2012[ ): defining V max as the 
maximum circular velocity of a halo and r max as the radius 
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Figure 4. Top panels: median formation redshift-mass relation for haloes identified at zq = in the ACDM simulation. The shaded 
region encloses the quartiles of the distribution at fixed halo mass. The solid line represents the prediction derived using the modified 
version of the |Giocoli et H] ( |2012b| > mass accretion history model. Bottom panels: differences with respect to the ACDM model for 
the three cDE models. The lines represent the predictions obtained from the modified model by |Giocoli et ah| ( |2012b| | by changing the 
definition of S(M) for each cosmology only. 



at which this velocity is attained so we can write: 



Ps 
Pcrit 



200 



c 200 



3 ln(l + c 2 oo) - caoo/(l + c 2oo) 

V s 2 

V max 

H(z) r m < 



= 14.426 



(19) 



where pcrit represents the critical density of the Universe. In 
Figure [7] we show the distribution function of the maximum 
circular velocity and of the radius at which it is attained for 
all haloes that at the present time are more massive than 
5 x 10 M(Tj/h. Since to faithfully determine the halo con- 



centration a minimum number of particles is required ( Neto 
et al.pOOT l, for this analysis we consider only haloes more 
massive than 5 x 10 13 Mq /h\ this lower limit ensures that 
the haloes are resolved with at least 850 particles. We notice 
that while the haloes in EXP003 and EXP008e3 have values 
of Vmax and Rmax that are not much different from the cor- 
responding haloes in the ACDM, the haloes in SUGRA003 



are more compact and typically have a higher maximum 
circular velocity that is attained at a smaller radius than in 
ACDM; as we will see later on, this corresponds to more con- 
centrated haloes in SUGRA003 at z — than in the other 
three cosmological models. 

In the context of the CDM hierarchical structure for- 
mation scenario, halo concentrations are thought to be rem- 
iniscent of the cosmic mean density at the time of collapse, 
thereby resulting in smaller objects having on average higher 
concentrations due to their earlier formation epoch. Consis- 
tently with this picture, clusters of galaxies have the lowest 
concentrations - typically of the order of ~ 4 - being the last 
collapsed structures in the Universe. At fixed redshift and 
halo mass the concentration tends to be larger for haloes 
that assemble their mass earlier; this also corresponds to 



haloes that are, on average, more relaxed. Neto et al. (20071, 
|Macci6 et aT] p007| [20(jt] > and |De Boni et al.| ( |2012[ ) have 
shown that including unrelaxed haloes in the sample results 
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Figure 5. Median mass accretion history of main halo progenitors. In the top part of each panel we show the results from the ACDM 
simulation for three different initial mass bins, where the shaded region encloses the quartiles of the distribution at fixed halo mass 
fraction assembled. In each panel, the solid curve represents the prediction from eq. JTsJ . The bottom part of each panel shows the 
residuals with respect to the ACDM simulation measurements of the three cDE models together with the residuals for the theoretical 
predictions changing S(M) in the mass accretion history model. Point style as in Figure H] 
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Figure 6. Residuals of the MAH measured in the three cDE models with respect to the MAH model changing the linear power 
spectrum normalization accordingly. The vertical dashed line in each panel shows the redshift at which the MAH model falls below the 
mass resolution equal to 5 X 10 12 Mq jh. 



in a c — M relation with a lower normalization and a larger 
scatter. 

Several studies based on the analysis of numerical sim- 
ulations have shown that halo structural properties are 
mainly related to the halo mass accretion history, and so to 
their environme nt ( |Gao et al.||2004| |Sheth fc Tormen|2004 
|Gao et aL"||2008[ ). Not only do less massive haloes possess a 
higher concentration but typically also a smaller mass frac- 



tion in substructures (De Lucia et al.||2004 van den Bosch 


|et al.||2005 


|Giocoli et al. 


|2010a|). In particular, the c — M 


models by 


Navarro et al. 


(1997|) and Bullock et al. (20011 



rely on the idea that the central density of the haloes reflects 
the mean density of the Universe at a time when the central 
region of the halo was accreting matter at a high rate. In 
this framework Zhao et al. ( 2009 1 define this epoch when 



the main halo progenitor is assembling 4% of its final mass. 

In order to understand how halo concentration corre- 
lates with the accretion history, in what follows we study 
the correlation between C200 and two conventional markers 



of the halo MAH: zo.5 i-e. the redshift at which the main halo 
progenitor assembles half of its final mass (we will denote 
with to. 5 the corresponding cosmic time) and Z0.04 (at which 
corresponds to. 04), i-e. the redshift at which it assembles 4% 
of its mass. 

In Figure [8] we show the correlation between the halo 
concentration and to /to. 5 where to represents the time at 
which the halo is considered - in our case zo = 0. In each 
panel we show the measurements for all haloes more mas- 
sive than 5 x 10 13 M Q //i in the ACDM and the three cDE 
models. The data points with error bars show the median of 
the correlation and the quartiles of the distribution at fixed 
to /to. 5- In the top-left panel the solid black line shows the 
least-square fit to the ACDM measurements, and the dotted 
red curve represents the following best-fit power law-relation 
obtained by minimizing the scatter: 



C200 



1 + 



to 



1.5to.i 



l/S 



(20) 
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Figure 7. Distribution function of maximum circular velocity 
(top panel) and of the radius at which this is reached (bot- 
tom panel) for all haloes at z = with masses larger than 
5 X 10 13 Afo//i in the four cosmological simulations. 



Eq. [20] has the same functional form of the equation pro- 



posed by Zhao et al. (20091. These two relations are also 



overplotted in the other panels. In each panel the gray his- 
tograms on the x and the y axis represent the distributions 
of to/to.5 and C200, respectively. The different line types in 
the EXP003, EXP008e3 and SUGRA panels show the least- 
squares fit to the corresponding measurements. In Table [3] 
we summarize these results and give an estimate of the rms 
of the data points for the different numerical simulation re- 
sults defined as: 



\ 



E 



(log(ci 



lQg(c mo rfei)) 

N 



(21) 



Zhao et al. ( 2009 1 have shown that the halo concentra- 



tion has a strong correlation with the time to. 04 at which 
the main halo progenitor assembles 4% its mass, in the idea 
that typically haloes acquire a concentration of 4 when they 
assemble 4% of their mass. From their measurements, the 
best fit relation between the concentration and to. 04 reads 
as: 



to 



3.75to.o 



1/8 



(22) 



where c V i r represents the ratio between the scale radius r s 



and the virial radius R V i r . We recall that Zhao et al. (20091 
define R V i r such that it encloses an overdensity A v i r accord- 
ing to the spherical collapse model, with respect to which 
they also define the virial mass M v i r . In Figure [9] we show 
the median of the correlation between C200 and to/to. 04 for 
the measurements in the four N-body simulations used in 
this work. The short-dashed lines show equation (22 1, that 
we recall is valid for the c V i r and M v i r definitions. In the 
top-left panel the solid line shows our best relation obtained 



modifying the Zhao et al. ( 2009 1 fitting function to be valid 



for the C200 and the M200 definitions, which reads as: 



C200 



1 + 



to 



3.2to.o 



1/13 



(23) 



These two curves are also shown in the other three panels 
of Figure |9| The different curve types in each panel show 
equation ( |23[ ) renormalized (see Table [4| in order to best 
fit the data points. In Table [4] we summarize the results of 
Figure [9] together with the rms with respect to the different 
models. 

Using the MAH model previously presented in this work 
and considering the correlations between C200 and to. 5 and 
C200 and to. 04, it is possible to estimate the concentration- 
mass relation model, in the following way. Given a mass 
M200 at redshift zq, we can estimate the redshift at which the 
main halo progenitor assembles half (or 4%) of its mass by 



using equation ( 13 1, and then the corresponding cosmic time 



by integrating over the scale factor a — 1/(1 + z) the inverse 
of the Hubble constant 1/H(a), from which we can obtain 
C200- In Figure, [To] we show the median of the concentration- 
mass relation for all haloes more massive than 5 x 10 13 Mq /h 
at redshift zo = 0. The data points in both right and left 
panels are the same. In the top panels they represent the 
median of the concentration-mass relation measured in the 
ACDM simulation, and the shaded region encloses the first 
and the third quartiles of the distribution at fixed mass. In 
the bottom panels we show the differences of the median 
of the measurements in the three cDE models with respect 
to the ACDM one. In Table [H we summarize the residuals 
of the c — M estimates with respect to the different model 
predictions. 



Maccio et al. (20081 andlZhao et al. (20091 have shown 



that haloes in standard ACDM simulations with low values 
of as and/or Q m tend to possess a lower concentration due 
to their later formation epoch. In what follows we will try 
to understand if this is also the case in our non-standard 
cDE simulations, i.e we will enquire whether the concentra- 
tion can be inferred from the ACDM one by simply taking 
into account the different normalization of the linear matter 
power spectrum at z = within the standard MAH model. 

From the analysis presented above, we have already no- 
ticed that even if the EXP003 and EXP008e3 models have 
a higher value of as, the halo concentration is not too dif- 
ferent from that measured for haloes of the same mass in 
the ACDM model. Only the SUGRA003 haloes are found 
to have a significantly higher concentration, irrespectively 
of their value of as that is very similar to the one in the 
ACDM run. These results are consistent with the previous 



findings of Cui et al. (20121 



In order to quantify to what extent halo concentrations 
are expected to change for different values of as within a 



Mass accretion history in cDE models 11 



Table 3. Rms of the data with respect to the models of Figure [8] 



model rms cq. ( 20 1 rms ACDM least-squares rms least-squares 



least-square values 



ACDM 
EXP003 
EXP008e3 
SUGRA003 



0.128 
0.132 
0.133 
0.255 



0.127 
0.134 
0.135 
0.254 



slope: 0.870, zero point: 0.503 
0.128 0.748, 0.506 

0.128 0.749, 0.488 

0.135 0.897, 0.704 




Figure 8. Correlation between the concentration and the formation time when the main halo progenitor assembles 50% of the mass at 
20 = 0. The four panels refer to the different cosmological models. In each panel the points with the error bars show the median with 
the two quartiles. Solid (black) and dotted lines in all panels are the same and represent the least-square fit to the haloes in the ACDM 
simulation and equation l |20| respectively. The different line types in the panels referring to the cDE models represent the corresponding 
least-square fit to their data points. 



ACDM cosmology, we plot in Figure[TT]the concentration-crs 
prediction combining the least-squares fit C200 — to/to.5 and 
the MAH model for a halo with mass M200 = 1 x 10 14 M Q //i 
as expected for ACDM - with the same cosmological pa- 
rameters of our ACDM run - but with different values of 
as- We take into account the change of as by renormalizing 



the mass variance such that: 

S CT8 (M) = 5 J|^Sacdm(M). (24) 

Fixing all the other cosmological parameters, we notice that 
the larger <rg the larger is the halo concentration, as ex- 
pected, reflecting the fact that haloes tend to form earlier. 
From the figure we can see that when as changes from 0.75 
to 1 the halo concentration changes by ~ 20%. The data 
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Table 4. Rms and least-squares fit parameters of the data with respect to the models shown in Figure [9] 



model 


Zhao et al. 


( 


2009 


\ eq. ( 


23 


1 cq. ( 


23 


renorm. 


re-normalization in log space 


ACDM 
EXP003 
EXP008e3 
SUGRA003 


0.152 
0.178 
0.168 
0.211 




0.143 
0.148 
0.147 
0.241 




0.143 
0.143 
0.143 


-0.04 
-0.05 
0.19 



Table 5. Rms of the data with respect to the models in the concentration-mass relation, as presented in the panels of Figure [Tp] 



model 


from eq. 


20 


I from 1-s c 2 oo-V*o.5, ACDM 


from 1-s C200-W*0.S 


from eq. ( 


23 


from eq. ( 


23 


1 renorm. 


ACDM 
EXP003 
EXP008e3 
SUGRA003 


0.153 
0.157 
0.156 
0.270 




0.154 
0.161 
0.161 
0.259 


0.153 
0.152 
0.162 


0.154 
0.165 
0.163 
0.260 




0.154 
0.151 
0.162 
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Figure 10. Concentration-mass relation at 20 = 0. The data points show the median c — M relation measured in the ACDM simulation 
and the shaded region encloses the quartiles. In the bottom part of each panel the points show the differences of the c — M measured in 
the three cDE simulations with respect to the ACDM one. In the left-top panel the two solid curves show the c — M prediction using 
equation \20\ and the least-squares relation C200 — to/*0.5- In the bottom left panel the curves show for each model the residuals of the 
c — M prediction from the least-squares fit in C200 — *oAo.5 relation with respect to the one for ACDM. In the top-right panel the c — M 
prediction is from equation ||23| and in the bottom panel the curves are the residuals of the predictions from renormalizing eq. |23|. 




0.75 0.8 0.85 0.9 0.95 1 

Figure 11. Concentration-erg relation for a halo with mass 
M200 = 1 x 10 14 Mq//i using the least-squares fit C200 — W*0.5- 
The colored data points show the median of the measurements in 
the simulations with the quartiles. For each simulation, we con- 
sider haloes within a logarithmic mass bin d log(M) = 0.001 cen- 
tered in M200 = 1O 14 M0 /h. The three black data points labelled 
as M-WMAP1, M-WMAP3 and M-WMAP5 are the estimates 
obtained using the best fit relation at 2 = by |Macci6 et al.| 
( 2008 ) for the three corresponding numerical simulations. 
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Figure 12. Subhalo mass function measured in the four cosmo- 
logical simulations used in this work. The subhalo mass is rescaled 
with respect to the host halo mass M20O1 an d for each halo we 
have considered only subhaloes with a distance from the host halo 
centre smaller than i?200- I n the bottom part of the figure we show 
the differences of the subhalo mass function with respect to the 
measurements done in the ACDM simulation. 

points with the error bars show the median and the quar- 
tiles of C200 for the same mass haloes measured in the four 
simulations at z = 0; we consider a logarithmic mass bin 
dlog(M) = 0.001 centered in Af 20 o = W 14 M e /h. The three 
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Figure 13. Mass fraction in substructures within R200 as a 
function of the host halo mass. In the bottom part we show the 
differences of the mass fraction in substructures with respect to 
the ACDM measurements. 



data points M-WMAP1, M-WMAP3 and M-WM AP5 are 
the predictions based on the C200 — A/200 relations by |Maccio| 
et al. ( 2008 1 for the three cosmological models at z — 0. 



The 

fact that the predictions from the relations by |Maccio et al. | 
( 2008 ) shift down with respect to the solid curve reflects the 
smaller values for Q m adopted their simulations: f2 m ,ACDM > 
Qm,M-wMAi > ft>m,M—wMA5 > flm,M-wMA3- From the 
figure we notice that haloes in EXP003 possess a higher 
concentration than in ACDM, consistently with the higher 
value of os(zq) that characterizes the model (see Table [TJ. 

However, the expected trend of C200 as a function of 
cr g is not followed both in EXP008e3 and SUGRA003. This 
point represents a central result of the present work, and 
deserves a more detailed discussion. It is well known that 
standard cDE models are characterized by a suppression 
of the non-linear matter power spectrum with respect to 
what ACDM predicts for any given value of as- This has 
been demonstrated both by running cDE N-body simula- 
tions with the same erg (as done e.g. in Baldi|2011a I and by 
comparing the non-linear matter power spectrum extracted 
from cDE N-body simulations normalized at zcmb ~ 1100 



to the predictions of HALOFIT (Smith et al. 20031 for the 



corresponding values of erg (see e.g. Baldi 2012b I. Such sup 



pression is determined by the action of the "friction term" 
(see the discussion in Section[3]above) that for standard cDE 
models has the effect of suppressing the growth of nonlin- 
ear structures. Correspondingly, for a fixed value of crg CDM 
haloes are found to be less concentrated in a cDE cosmology 



when compared with respect to ACDM (see e.g. Baldi et al. 
|2010a| |Baldi||2011a| |Li fc Barro^[20TT| . For more complex 
cDE models, like time-dependent couplings or the "bounc- 
ing" cDE scenario, the role played by the "friction term" is 
not so straightforward as the latter can change sign during 
the cosmic evolution or have a non-trivial interplay with the 



time evolution of the coupling itself. In such more general 
scenarios, halo concentrations at fixed crg can be both higher 
and lower than in ACDM, depending on the specific model. 

Among the models considered in the present work, 
the slight increase of halo concentrations for the standard 
EXP003 cDE model shows how the suppression associated 
to the "friction term" is compensated by the higher value 
of crg of this cosmological model as compared to ACDM. 
Furthermore, the lower and higher values of halo concen- 
trations with respect to what is predicted by the standard 
ACDM scenario, that are observed for the EXP008e3 and 
SUGRA003 models, respectively, show how the dynamical 
evolution of the DE scalar field can alter the structural prop- 
erties of CDM haloes in a way that is clearly independent of 
the evolution of linear density perturbations. For instance, 
in the case of the "bouncing" cDE model (SUGRA003), halo 
concentrations significantly grow in time at low redshifts (as 
already found by Cui et aL]|2012 I due to the particular dy- 
namics of the DE scalar field that inverts its direction of 
motion at Z[ nv ~ 6.8 the reby also cha nging the sign of the 
"friction" term fi c 4> (see 



Baldi 



2012a 



for a more detailed 
discussion of the dynamics of "bouncing" cDE). The latter 
then acts as a dissipation term for virialized objects inducing 
an adiabatic contraction of halos that consequently evolve 
towards more concentrated virial configurations. 

This direct relation between the dynamics of the under- 
lying DE field and the formation and evolution of nonlinear 
structures offers the interesting prospect of using the forma- 
tion history of CDM haloes and their structural properties 
to disentangle cDE effects from possible variations of cos- 
mological parameters (as e.g. crg), of the linear galaxy bias, 



and of the mass of cosmic neutrinos (see e.g. La Vacca et al. 
2009} |Marulli et al.|2011 |. 



6.2 Halo substructures 

Studying the subhalo population in DM haloes extracted 
from a standard cosmological ACDM simulation at z = 0, 



Giocoli et al. (2010a I found that more concentrated haloes, 



forming at higher redshift, tend to possess on average less 
substructures above the same minimum mass. They have 
also stressed the fact that considering haloes with the same 
mass, the ones that form earlier not only possess a higher 
concentration, but also less subhaloes. This is because halo 
progenitors are accreted earlier, then spending more time 
in the potential well of the host halo, and thereby tend to 



lose a larger fraction of their initial mass (van den Bosch 
et al.|[2005 Giocoli et al.|[2008 |. To see if this phenomenol- 
ogy is also reflected in the cDE simulations considered in 
this work, we plot in Figure [12] the subhalo mass function 
considering for each halo all substructures resolved by sub- 
find with a distance from the host halo centre smaller than 
^200- In the top part of the figure we show the subhalo 
mass function for three different host halo mass bin (as in 
Figure [5| at z = 0, while in the bottom one the differences 
with respect to the measurements in the ACDM simulation. 
Since we have shown that haloes in SUGRA assemble ear- 
lier and are more concentrated than those in the ACDM 
run, we would have expected them to possess less substruc- 
tures, on the contrary to what is shown in the figure. Also 
the subhalo population in the EXP003 and EXP008e3 cos- 
mologies does not reflect their mass accretion history. In 
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these models, haloes possess less substructures than those 
in ACDM since they form earlier, but even if in EXP003 
haloes form earlier than those in EXP008e3 they are still 
found to host more substructures than in the latter model. 
This statement is further demonstrated in Figure |l3| where 
we show the halo mass fraction in subhaloes as a function 
of Af2oo- In the bottom part of the figure we show again the 
differences of the measurements in the cDE with respecto to 
the one in ACDM: haloes in SUGRA (EXP008e3) are more 
(less) substructured than those in ACDM by ~ 15% while 
the difference between EXP003 and ACDM is only of the 
order of a few percent. We argue that the main reason why 
SUGRA003 (EXP008e3) possesses more (less) substructures 
at z — is that halo progenitors at any redshift z < 1, which 
end up in present-day substructures, are more (less) concen- 
trated than in ACDM; so the gravitational heating and tidal 
stripping are less (more) efficient in disrupting the satellites. 



7 SUMMARY AND CONCLUSIONS 

Using state-of-the-art cosmological simulations for cDE 
models from the CoDECS Project we have studied how col- 
lapsed haloes grow as a function of the cosmic time, and how 
present-day systems acquire their structural properties. We 
summarize our study and main results as follows. 



We updated the model developed by |Giocoli et al 



( 2012b I in the context of the extended- Press & Schechter 
( 1974[ ) formalism to be able to describe the halo mass ac- 
cretion history of a ACDM cosmology when haloes are iden- 
tified by their M200 mass (instead of M v i r ), defined as the 
mass of the spherical region around the halo center enclosing 
an average density 200 times larger than the critical density 
of the Universe. 

• By rescaling the MAH model with the normalization 
of the linear matter power spectrum for the different cDE 
runs, we noticed that the simulation results are reproduced 
quite well for small masses and low redshift. For galaxy 
cluster-size haloes and up to high redshifts, when less than 
1% of the present-day mass is assembled, the models and 
the simulation measurements differ by about 20%. Haloes 
at z = in the different cDE cosmologies (standard cDE, 
time-dependent cDE, and "bouncing" cDE) typically assem- 
ble their mass at higher redshift as compared to those in the 
standard ACDM run. 

• We studied the formation redshift zj for different frac- 
tions / of the assembled mass: standard and time-dependent 
cDE show a systematically higher Zf by ~ 20 and ~ 10%, re- 
spectively, as compared to ACDM. On the other hand haloes 
in our "bouncing" cDE cosmology have a formation redshift 
quite similar to those in ACDM, for small /, while for large 
assembled fractions Zf can also be higher than in ACDM by 
up to ~ 20%. 

• Analyzing the correlation between the halo concentra- 
tion and two typical formation times to. 5 and to. 04, we con- 
firm that the usual correlation between halo formation time 
and halo concentration at z — (the earlier the formation 
time, the higher the concentration) still holds in cDE cos- 
mologies. Generalizing the prediction for ACDM we provide 
some fitting functions for such correlations also in cDE mod- 
els. In particular, in our "bouncing" cDE scenario haloes are 



very concentrated at z = 0, inconsistently with the corre- 
lation expected for standard cosmologies ( Cui et aL~||2012 1. 
Such structural properties could produce very compact and 
luminous galaxies located in the centre of haloes and also 
galaxy clusters that are very efficient for strong lensing. 

• Considering the correlation between C200 and the time 
at which the main halo progenitor assembles half of its mass, 
we have confirmed that for standard ACDM cosmologies the 
concentration is a monotonic function of as'- in particular for 
a cluster-size halo the change in concentration between two 
ACDM cosmologies with ag — 0.75 and as = 1 is of about 
20%. Moving to cDE cosmologies, we have shown that only 
the standard cDE run, characterized by an exponential self- 
interaction potential and a constant coupling, is found to be 
in agreement with such predictions, while both cDE models 
with a time-dependent coupling and the "bouncing" cDE 
model are not found to be consistent with the expected evo- 
lution of halo concentrations with as- This inconsistency 
offers the interesting prospect of disentangling the effects of 
dark interactions from the variation of standard cosmolog- 
ical parameters using the evolution of halo concentrations 
HGiocoli et"aL|2012c| ). 

• The standard mass accretion history model is also found 
to be not directly applicable to predict the subhalo popu- 
lation of the host haloes at z = within cDE cosmologies. 
In fact, if we generalize to the cDE simulations the state- 



ments by Gao et al. (20041, van den Bosch et al. (20051, 



and Giocoli et al. ( 2008 1 - valid for standard ACDM - that 



haloes with a higher formation redshift typically possess less 
substructures, we would expect to find more substructures 
in ACDM than in all our cDE models at any fixed halo 
mass. On the contrary, we find that the trend is actually 
reversed, with haloes of the same mass having ~ 15% more 
substructures within -R200 in the "bouncing" cDE run than 
in ACDM, while haloes forming within a time-dependent 
cDE cosmology have less substructures by the same amount, 
while haloes in standard cDE are roughly consistent with 
their higher formation redshift and concentration but have 
slightly less substructures then in ACDM. We argue that 
the higher concentration and the larger number of clumps 
in the SUGRA003 model are due to the fact haloes formed 
when the average density of the universe was larger, and the 
difference in the average density between the formation time 
of the small halos and the larger ones was also larger than in 
the in ACDM. This makes satellites not only more concen- 
trated than in thoes in ACDM, but also more concentrated 
relative to the host halos: the subhalos then tend to survive 
tidal striping for longer time. In summary, while in stan- 
dard ACDM simulation - indepdentently of the small scale 
behavior of the linear power spectrum used to generate the 
initial conditions ( Schneider et al. 2012 1 - the concentration- 



mass relation is expected to prove the halo formation time, 
the subhalo abundance as a function of the host halo mass 



validate the dynamical friction and the tidal stripping (van 
|den Bosch et al.|2005| ). 



To conclude, in the present work we have performed a 
detailed analysis of the mass accretion history of collapsed 
haloes in a sample of coupled Dark Energy cosmologies in- 
cluding different choices of the Dark Energy self-interaction 
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and coupling functions. We have studied how haloes acquire 
their structural properties along the cosmic time and tested 
whether it is possible to attribute the detected differences 
with respect to the standard ACDM case to the different lin- 
ear power spectrum normalization. Interestingly, we found 
that this is not possible for all coupled Dark Energy models, 
and we identified which observables allow to break such a 
degeneracy. Finally, we have investigated by how much halo 
concentration and subhalo abundance deviate from ACDM 
for the three coupled Dark Energy models included in our 
simulations set. In particular, we showed that the unex- 
pected high concentration and dumpiness of haloes make 
the "bouncing" coupled Dark Energy model particularly in- 
teresting in the context of both weak and strong gravita- 
tional lensing. We intend to study the possibility of distin- 
guishing such models from ACDM using lensing data in fu- 
ture work. 



APPENDIX A: MASS ACCRETION HISTORY 
IN TERMS OF THE VIRIAL MASS FOR THE 
ACDM SIMULATION 

In this Appendix we show that the MAH model built follow- 
ing back in time the haloes in the ACDM simulation in terms 
of their virial mass M„; r is in perfect agreement with the 

we 



results obtained by Giocoli et al. ( 2012b I. In Figure Al 



show the cumulative generalized formation redshift distribu- 
tion when 90%, 50%, 10% and 1% of the halo mass is assem- 
bled in term of the universal variable Wf. We compute the 
redshift at which the main halo progenitor assembles a frac- 
tion / of its mass at z = interpolating its mass accretion 
history along the different simulation snapshots, and then 
compute Wf = (5 c (zf) — 8 c (z )) /y/S{fM vi . 



The solid line in the figure shows equation ( 12 1 with the 



S(M vi 



best fit value of a/. The dashed curve refers to the relation 
by Nusser & Sheth ( 1999 1 modified as proposed by Giocoli 



et al. (20071 and the dot-dashed curve equation ( 12 I with 



the cxf — f relation by Giocoli et al. ( 2012b I. We notice that 
our best fit is in very good agreement with what found by 



Giocoli et al. (2012b I, and does not depend on the differ- 



ent cosmological parameters of the simulation since Wf is 
a universal variable (Press & Schechter 1974 Bond et al. 
1991| |Lacey fc Cole||1993[ |1994[ ). The small differences can 
be traced back to the different code used to run the two 
numerical simulations and to do the post processing analy- 
ses. While HYDRA (|Couchman et aLfll995l ) and GADGET 



(Springel et al. 2001a I have been used for the GIF2 sim- 
ulation studied by Giocoli et al. (2012b I, a modified ver- 



sion of GADGET2 ( |Sprin gcl 2005) developed to include all 



the additional physical effects that characterize cDE models 



(Baldi et al. 2010b I has been run four our ACDM simula- 
tion. For the post-processing analyses Giocoli et al.| ||2012b l 



have adapted and run the pipeline presented by |Tormen 
et al.| (20041, while here we have run the codes by Springel 



et al. 



(2001a I and Boylan-Kolchin et al. (20091. We recall 



the reader also on the fact that the GIF2 is a pure DM sim- 
ulation while the CoDECS presents a baryon fluid with no 
hydrodynamic treatment included. 

In Figure |A2| we show the correlation between the free 
parameter of the MAH model aj and the assembled mass 
fraction. The data points represent the best fit values ob- 




0.4 0.6 
fraction 

Figure A2. As Figure[3]but in terms of M v i r . 

tained fitting the cumulative formation redshift distribution 
for different values of /; the shaded region encloses la and 
2a contours around them. In the figure, the solid curve rep- 
resents the following relation: 

0.837 



a f = 



f° 



while the dot-dashed one is for: 
0.815 



a f = 



-2f 



-2f<> 



0.707 



(Al) 



(A2) 



as found by Giocoli et al. ( 2012b I. The two relations for the 



ACDM runs are in perfect agreement confirming the fact 
that the MAH model can be generalized for any ACDM 
cosmology, independently of the cosmological parameters. 



APPENDIX B: PUBLICLY AVAILABLE 
MERGER TREE FILES ON THE CODECS 
DATABASE 

The merger trees of the different L-CoDECS simulations used 
to perform the analysis discussed in the present paper have 
been produced using a linking algorithm outlined in |Spring"eI| 
|et al.| ( |2005[ ) and |Springel et al.| ( [2008| > . With such algorithm, 
we produced a single merger tree file for each cosmological 
model of the L-CoDECS suite (i.e. even for those models that 
have not been considered in the present paper). 

As an extension of the public CoDECS database, 
we hereby release the merger tree files that are 
now directly available at the CoDECS website 
( [http://www.m arcobaldi.it/CoDECS). These are un- 
formatted binary files with an average size of about 10 Gb, 
and detailed instructions on how to read and use the data 
can be found on the CoDECS guide (version 2.0) that can 
also be directly downloaded from the CoDECS website. 

The access to these files is subject to the same terms of 
use that apply to the whole CoDECS public database. 
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Figure Al. As Figure [2] but in terms of M v i r 
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